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Abstract 

A statistical study is carried out on the photospheric magnetic nonpotentiality 
in solar active regions and its relationship with associated flares. We select 2173 
photospheric vector magnetograms from 1106 active regions observed by the 
Solar Magnetic Field Telescope at Huairou Solar Observing Station, National 
Astronomical Observatories of China, in the period of 1988-2008, which covers 
most of the 22nd and 23rd solar cycles. We have computed the mean planar 
magnetic shear angle (A</>), mean shear angle of the vector magnetic field (A-0), 
mean absolute vertical current density (| J z \), mean absolute current helicity den- 
sity (|/i c |), absolute twist parameter (|a a v|), mean free magnetic energy density 
(Pfrco), effective distance of the longitudinal magnetic field (c?e), and modified 
effective distance (<iEm) of each photospheric vector magnetogram. Parameters 
|^c|, Pfree, and d,Em show higher correlation with the evolution of the solar cycle. 
The Pearson linear correlation coefficients between these three parameters and 
the yearly mean sunspot number are all larger than 0.59. Parameters A(f>, Ai/j 7 
\Jz\, |cy av |, and g?e show only weak correlations with the solar cycle, though the 
nonpotentiality and the complexity of active regions are greater in the activity 
maximum periods than in the minimum periods. All of the eight parameters show 
positive correlations with the flare productivity of active regions, and the com- 
bination of different nonpotentiality parameters may be effective in predicting 
the flaring probability of active regions. 
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1. Introduction 



It is generally accepted that the magnetic field dominates the evolution of the 
Sun, and the magnetic nonpotentiality may be responsible for massive energy 
released in solar explosive activities such as flares and coronal mass ejections 
(CMEs). Even though the physical parameters that characterize the nonpoten- 
tiality of the magnetic fields do not directly trigger solar eruption, active regions 
(ARs) with strong nonpotentiality and considerable complexity are more likely 
to explode than simpler ones because the former have sufficient free energy to 
release. Studying the evolution of magnetic nonpotentiality and its relationship 
with solar eruption from a statistical point of view is useful in predicting eruptive 
events and in understanding the long-term evolution of solar activity. 

The characteristics of nonpotential magnetic fields associated with flares in 
ARs have been studied for a couple of decades. Several nonpotential parameters 



of flare-productive ARs are investigated by many authors in detail. Gary et 



(1987) analyzed nonpotential features based on the observed vector magnetic 
field of NOAA AR 2684, and tried to answer the questions such as the relation 
between the nonpotential features and the flares. Wang et al. (|1996|) examined in 
detail the relationship between flare occurrence and nonpotentiality development 
in a flare-productive region. IMoon et al.\ (|2000|) analyzed the spatial distribu- 
tion and temporal evolution of several nonpotential parameters in NOAA AR 
5747. [Deng et ai.|(|2001) studied the daily evolution of magnetic nonpotentiality 
focusing on NOAA AR 9077 that produced a giant flare from 11 to 15 July 
2000. ILeka and B arnes (2003]) investigated the magnitudes and temporal varia- 
tions of several photospheric magnetic parameters, such as horizontal magnetic 
gradient, vertical current, current helicity, twist parameter a, magnetic shear 
angles, and excess magnetic energy in three ARs, to identify properties that are 
important to the solar flare production. IDun et ai\ (|2007[) calculated the mag- 
netic shear angle, vertical current, and current helicity of NOAA AR 10486 and 
analyzed their temporal evolutions and spatial relations to the flares. Falconer, 
Moore, and Gary p002l 12006) also discussed the relationship between several 
nonpotcntiality-rclated parameters and CME productivity of ARs. Some authors 
studied the physical characteristics of photospheric magnetic fields in ARs and 
their relationship to flares in a statistical sense. ICui et ai\ (|2006[) preferred the 
maximum horizontal gradient, length of the neutral lines, and number of singular 
points; |Jing et ai\ (|2006) used the mean value of magnetic gradients at strong- 
gradient magnetic neutral lines, length of strong-gradient magnetic neutral lines, 
and total magnetic energy; ICui et al\ (|2007|) also tried the length of neutral lines 
with strong gradient and/or strong shear; [Cui and Wang (2008) studied the total 
unsigned current and total unsigned current helicity. IGuo et aT\ (|2006l 120071 
12010) utilized the effective distance that could roughly quantify the magnetic 
complexity of ARs, to study the evolution of magnetic complexity in ARs, its 
relationship with erupting activity of ARs, and its evolution with the solar cy- 
cles, respectively. A close relationship between magnetic nonpotentiality of ARs 
and solar eruptions is observationally studied in numerous articles. IRust et al.\ 
(|1994) reviewed the progress of pre-flare state studies, including nonpotential 
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characteristics, and pointed out the importance of the more quantitative results 
from the vector magnctographs. 

Lacking of long-term steady vector magnetic field observations, the line-of- 
sight magnetograms are mostly used as means of case analyses and statistical 
researches. It is necessary to study the long-term evolution of appropriate param- 
eters derived from vector magnetic fields and find whether there is a correlation 
between these parameters and flares or other eruption events. So far ground- 
based telescopes, whose maintenance for long-term observations is relatively easy, 
still have the advantages in data accumulation. Since 1987, the Solar Magnetic 
Field Telescope (SMFT) QAi and Hu, 1986D , installed at Huairou Solar Observing 
Station (HSOS), National Astronomical Observatories of China (NAOC), has 
reliably worked for more than twenty years that almost covers two complete 
solar cycles. Based on the long-term observations of SMFT, a statistical study 
on photospheric magnetic nonpotential parameters in ARs is carried out, and 
the relationship between these parameters and solar flares is also presented. 

In Section [3J observations by SMFT and data reductions are introduced. 
Measures of nonpotcntiality and complexity arc shown in Section [31 In Section 
21 statistical analysis and results on magnetic nonpotentiality associated with 
flares are presented. Conclusions are given in Section [5] 

2. Observations and Data Reductions 

2.1. Observations at Huairou Solar Observing Station 

SMFT is a 35cm-apcrturc filter-type vector magnctograph that measures solar 
magnetic and velocity fields in both the photosphere and chromosphere. The 
magnetograms used in this study are obtained by one of its working spectral 
lines, Fe I A5324.19 A that measures photospheric vector magnetic fields. The 
Fe I A5324.19 A line is a strong and broad line with an equivalent width of 
abo ut 0.334 A and a Lande factor of 1.5 (|Ai, Li, and Zhang, 1982| Wang et ai, 
1996). Measuring at an offset -0.075 A fr om the line center of Fe I A5324.19 
A, SMFT der ives the longitudinal component of the vector magnetic field from 
the circular polarization (Stokes parameter V), and the transverse components 
from the linear polarization (Stokes parameters Q and U) at the line center. It 
takes about 3 min to get the data for generating a vector magnctogram. The 
observations aim at limited areas of the solar disk, being mainly the ARs where 
sunspots are located. From 1988 to 2008, three different CCDs were used in 
SMFT, and the data were captured by them with three different sizes of the 
field of view (FOV): before 25 August 2001, the FOV was 5.23' x 3.63' and the 
pixel scale was 0.61" x 0.43"; from 25 August 2001 to 30 November 2001, the FOV 
was 4.06' x 2.77' and the pixel size was 0.48" x 0.32"; after 1 December 2001, the 
FOV was 3.75' x 2.81' with the pixel resolution of 0.35" x 0.35". To improve the 
signal-to- noise ratio (SNR), the observing system adopts a 256- frame integration 
for routine observations. After smoothing the three sets of data with 3x5,3x5, 
and 5x5 pixels, respectively, we get a better SNR and the spatial resolution has 
been reduced to about 2" x 2". The 3er noise of a vector magnetogram is about 
20 gauss (G) for the longitudinal field and 150 G for the transverse field. 
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Table 1. Number distributions of selected magnetograms and active regions from 
1988 to 2008. The counts of the magnetograms (MGs) and active regions (ARs) in 
the solar maximum periods are relatively higher than those in the solar minimum 
periods during the same solar cycle. The samples have reasonably covered two solar 
cycles. 





1988 


1989 


1990 


1991 


1992 


1993 


1994 


1995 


1996 


1997 


1998 


MGs 
ARs 


10 

8 


37 
21 


47 

30 


108 
58 


138 
65 


97 
52 


74 
39 


56 
25 


11 

6 


54 
32 


114 
62 






1999 


2000 


2001 


2002 


2003 


2004 


2005 


2006 


2007 


2008 


Total 


MGs 
ARs 


166 
88 


361 
166 


285 
136 


233 
123 


144 
77 


67 
34 


101 

50 


46 
22 


20 
11 


1 

2 


2173 
1106 



The vector magnetograms obtained from various magnetographs of different 
observatories were compared for the research of magnetic field structures in 
specific ARs ( |Wang et al, 1992[ |Bao et al, 2000} |Zhang et al, 2003D and for 
the statistical analysis ( |Pevtsov, Dun, and Zhang, 2006} |Xu et al, 2007[ ), to 
confirm the accuracy of the magnetic field measurement by SMFT. As a result, 
these data are basically consistent, including the longitudinal and transverse 
components of the magnetic field, electric current, helicity parameters, and so 
on. Small discrepancies due to observing methods and data inversion methods 
are insignificant in the present statistical study of magnetic nonpotcntiality. 



2.2. Data Selection and Preprocessing 



The samples located within 30 from the solar disk center have been selected to 
minimize the influence of projection effects. Of all the vector magnetograms, only 
one magnetogram is picked up for an AR each day. Most of selected samples are 
observed during 01:00 to 07:00 UT when HSOS is in a relatively better seeing 
condition. 2173 observed photospheric vector magnetograms containing 1106 
ARs from June 1988 to March 2008 are chosen as samples. Table [T] lists the 
selected magnetograms and ARs each year in the two cycles. Various typos of 
ARs arc included in this sample table, from simple unipolar regions to complex 
5 ARs. From 1988 to 2008, there are 6095 ARs recorded by NOAA, and the 
percentages of the ARs that produced >C-, >M-, and X-class flares are 39.1%, 
15.0%, and 1.0%, respectively. Our selected ARs account for 18.1% of the total 
ARs over this period, and 71.3%, 20.7%, and 1.0% of the selected ARs produced 
>C-, >M-, and X-class flares. Though these ARs are biased to the larger and 
more active regions, they only produced 18.6% of all the AR-related >C-class 
flares (C-class flares, 18.8%; M-class flares, 17.2%; X-class flares, 19.8%) within 
the subsequent 48 h from the time of the selected magnetogram. Therefore, the 
sample coverage is reasonable and representative. 

The one-to-one match between Huairou and NOAA AR numbers is manually 
done for analyzing the relationship between the magnetic characteristics of an 
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AR and associated flares. Thus the flare records taken by Geostationary Opera- 
tional Environmental Satellites (GOES) can be used in the study. As the FOV 
of SMFT is fixed, each magnetogram containing a single complete AR near the 
center of the FOV is sampled. Therefore, excluding the possible appearance 
of other ARs in a particular magnetogram, the data reduction is uniformly 
carried out on all the selected magnetograms. The lists of soft X-ray (SXR) 
flare events are recorded by GOES in detail, which can be downloaded from 
National Geophysical Data Center (NGDCfl- 

A vector magnetogram consists of the longitudinal (Bi) and transverse (B t ) 
components of the magnetic field, after processing the raw data of Stokes /, Q, 
U, and V signals. The calibration method commonly used for SMFT is based 
on the relationship 

Bt = C t j, B t = C t {j[^) = (1) 

where /, Q, U, and V are Stokes parameters, and Cg, and Ct are the correspond- 
ing calibration coefficients of Bg and B t . The azimuth angle (</>) of the transverse 
field is obtained by 

= \ arctan (|) . (2) 

Fro m several calibrations on SMFT vector magnetographs ( Ai, Li, and Zhang, 
1982;|Wang, Ai, and Deng, 1996||Su and Zhang, 2004|), C t = 8381 and C t = 6790 




( |Su and Zhang, 2004 ) are selected in the present study. The statistical results 
are not affected significantly by the coefficient selection. 

The potential fie ld approximation method (|Harvey, 1969 Sakurai, Makita, 



and Shibasaki, 1985) is applied to resolve the 180 ambiguity of the transverse 
field. This method works well for the regions where the actual magnetic shear 
angles are less than 90 . Though the azimuth ambiguity resolution would be 
rough in a few regions that arc very nonpotential, the analyzed parameters 
we use are all macroscopic and averaged quantities over a whole AR, thus the 
influence could be neglected in the statistical analysis. All the magnetograms 
are not deprojected because the selected samples are located near the solar disk 
center ( |Gary and Hagyard, 1990] |Li, 2002||Cui et ai, 20071 ). 



3. Magnetic Nonpotentiality and Complexity Parameters 

3.1. Magnetic Shear Angle 

The magnetic shear is one of the typical parameters to describe the nonpoten- 



tiality of magnetic fields in the solar photosphere. |Hagyard et aT] (|1984p and Lii, 



Wang, and Wang (jT993p introduced the planar and three-dimensional magnetic 



^ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SOLAR_FLARES/XRAY_FLARES/| 
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shear angles, respectively. The three-dimensional magnetic shear angle Atp is 
the angle between the directions of the observed vector magnetic field B a and 
its corresponding potential field B p . The planar magnetic shear A<j), which can 
be considered as the projection of Atp onto the plane of the sky, is defined as 
the azimuthal difference between the observed field and the potential field. Atp 
is a more direct measure of deviation from a potential field than A0, and has 
a direct relation to the magnetic free energy of the nonpotential fields as well 
(see Section [331 |Lii, Wang, and Wa ng, 1993]). Here Acj) and Aip, which are the 
angles between the directions of two vectors, are both without sign: 

A0 = (B3t P ) = arccos (j|r^|j) , (3) 



A^ = (B ,B p ) = arccos (]|rprf) , (4) 

where the subscripts o and p indicate the observed and potential magnetic fields, 
respectively, and t refers to the transverse component of the magnetic field. 

3.2. Vertical Current Density, Current Hclicity Density, and Twist Parameter a 

The electric current, current helicity, and twist parameter (force-free field factor) 
a are also important to describe magnetic nonpotcntiality. Early in 1960s, solar 
phys icists began to observationally study the electric currents in ARs ( Severny, 
1965) and their association with solar flares (|Moreton and Severny, 1968|. Wang, 
Xu, and Zhang Tl994[) and |Schrijver| (|2007p indicated that the emergence of 
mag netic fields carrying electric currents may cause most maior solar flares. B ao 
et aZ. (|1999|) testi fied that the current helicity fcK , Scchafc rlll990l Pevtsov. Can- 
ficld, and Metcalf TT995jT of ARs plavs an important role in solar flare evolution. 
|Nakagawa and Raadu| (|1972p theoretically described the force-free field factor 
a as an indicator of the twist of the magnetic field lines; they deduced that 
the magnetic energy increases with the value of a. INindos and An drews (2004) 
discussed CME-associated big flares of 133 events to explain the roles of the 
magnetic helicity and force-free factor a of the pre-flare corona. 

According to Ampere's law and neglecting the effect of the electric displace- 
ment current, the vertical component of electric current density is obtained 

by 

r 1 m v m 1 fdB y dB x \ 

■h = —(V x B) z = — - — — , 5 

Mo Mo \ ox dy J 

where J z is in units of A km~ 2 , /io = 47r x 10~ 6 G km A -1 is the permeability in 
free space, and B x and B y are the two perpendicular components of a horizontal 
magnetic field. 

The current helicity is defined as 



H c = [ B-{VxB) dV, 
Jv 
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and the vertical current helicity density can then be calculated by 

W,=»,(V» B ), = B ,(|-|). (6) 

So far there is no observational method to obtain the transverse component of 
current helicity density (h c )t, and the symbol h c is usually used to represent its 
vertical component. 

On the basis of the force-free field assumption, the Lorentz force is equal to 
zero, i.e. the electric current is parallel to the local magnetic field: 

(V x B) x B = 0, 



or 



VxB = olB. 



If the magnetic field is approximated to be a linear force-free field, the local twist 
a is expressed as 

(V x B) z 



B z 

We use a a v ( |Hagino and Sak urai, 2004[ ) to characterize the dominated twist 
status of an AR: 

_ E(V x B) z ■ sign[^] 
3.3. Free Magnetic Energy Density 

The energy released through solar flares and other explosive events relies on the 
accumulation of the free magnetic energy (nonpotential magnetic energy) that 
is defined as the difference between the total magnetic energy (E) and potential 
magnetic energy {Ep): 

AE = E - E p . 



Hagyard, Low, and Tandberg-Hanssen (|198ip introduced the source field to 



describe the nonpotcntiality of a. magnetic field on the photosphere: 

B s = B a — Bp, 

where B Q is the observed vector magnetic field, B p denotes the potential field 
extrapolated from the longitudinal component of B Q , and B s is the so-called 
source field which is the nonpotential component of the magnetic field. 
The magnetic energy density of the source field is proportional to B^ : 

B 2 S {B -Bpf 
Phcc ~ 8^ 8^ ' (8) 
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Equation (J5J is used to calculate the free magnetic energy density on the pho- 
tosphcric layer in the present work. It is deduced from Equation that 



Plr. 



(B Q - B p f 
8tt 



B B p . 2 
— - — - sin 

2tt 



¥) 



in which B Q = \B Q \ and B p = \ B p \ . We can then find that the shear angle of the 
vector magnetic field Aip is directly related to the nonpotential magnetic energy 



( jLii, Wang, and Wang, 1993] [Wang et al, 1996[ ). 
3.4. Effective Distance 



Effective distance (dp), a structural parameter of an AR. proposed bv Chu- 
mak and Chumak (j!987p . presents a distinction between flare-quiet and flare- 
imminent ARs (Chumak, Zhang, and Guo, 2004 1. As a quantified magnetic 



complexity, c?e depicts the degree of the isolation or mutual penetration of the 



two polarities of an AR in a geometrical sense QGuo et al., 2006[ ). is calculated 
from 



where 



Bp + Rn 
Rpn 



(9) 



A p (A n ) is the total area of positive (negative) polarity, R p (R n ) is the equivalent 
radius of positive (negative) polarity, respectively, and R pn is the distance be- 
tween the flux- weighted centers of the two opposite polarities (|Guo et all [2006 , 



[20071 120T01 [Guo and Zha"nj[ 120071) . 

Considering the limitation of the parameter g?e to uniquely characterize the 
complexity of an AR, an additional factor \B Z \ is multiplied to c?e- The modified 
effective distance d^m is defined as 



4m = dn\B z \ = 



Rp + Rn 
Rpn 



BA 



(10) 



g?e is dimensionless, but \B Z \ provides d^m some practical physical meaning to 
some extent. Here, the dimension of d^m is the same as that of the magnetic 
field, in the unit of gauss, d^m also reflects the complexity of an AR, and could 
be considered as the weighted complexity of entire AR in terms of the averaged 
strength of the sunspot magnetic field. 



4. Statistical Analysis and Results 

4.1. Strength Distribution of Nonpotentiality during Solar Cycles 22-23 

By means of the magnetic nonpotentiality and complexity parameters described 
in Section [21 Equations (|3^)- ([TU)) arc used to calculate, for each photospheric 



SOLA: S0LA1731vl.tex; 7 March 2013; 23:48; p. 8 



Statistical Study on Magnetic Nonpotcntiality and Associated Flares 




Figure 1. A vector magnetogram of NOAA AR 5356 at 02:12 UT on 15 February 1989. The 
white and black regions correspond to the positive and negative longitudinal magnetic fields, 
respectively. The arrows indicate the directions of transverse fields, the lengths of which denote 
their magnitudes. Only the arrow vectors with St > 50 G are shown in the figure. The white 
and black solid contours indicate B\ = ± 20 G; the white and black dashed contours indicate 
Bi = ±80 G. The red solid contours arc for B t = 200 G. The FOV is about 5.23' X 3.63'. 

vector magnetogram, the mean values of the planar magnetic shear angle (A(f>), 
shear angle of the vector magnetic field (Aip), absolute vertical current density 
( | J z | ) , absolute current helicity density ( | h c | ) , and free magnetic energy density 
(pfroo), and also calculate the absolute twist parameter (|a av |), effective distance 
of the longitudinal magnetic field (c?e), and modified effective distance (<iEm)- 
In each of the magnetograms, only the areas where the strength of longitudinal 
magnetic field is greater than 20 G arc used in the calculations for the parameters 
\J Z \, \h c \, |ct av |, pfroe, and d,Em- The two mean shear angles Atfi and Atp are 
obtained on the pixels with both the transverse magnetic field greater than 200 
G and the strength of longitudinal field greater than 20 G, which are distributed 
at the sunspot penumbra and along the polarity inversion lines. The longitudinal 
field strength of 80 G is chosen as the lower threshold for obtaining c?e- Figure 
[T] shows one of the magnetogram samples, which is NOAA AR 5356 observed at 
02:12 UT on 15 February 1989 by SMFT. The white and black solid contours 
(B, = ±20 G) depict the boundaries for calculating A<j>, Aip, \J Z \, \h c \, |a av |, 
Pfree, and 4m- The white and black dashed contours (B z = ± 80 G) enclose the 
regions for calculating d&. The red solid contours (B t = 200 G) mark the edges 
for further confining the pixels used to obtain A<f> and Aip. 

Figure [2] shows the distributions of the above eight parameters of ARs as 
a function of time from 1988 to 2008. For reference, two plots of the monthly 
mean sunspot numbers are displayed above other panels. The active samples 
(i.e., flare-productive ARs) arc defined as the ARs with the equivalent flare 
strength (i.e. flare index FI defined in Section [472]) greater than a typical value 
(such as M1.0 here) within the same subsequent time window (24 h in Figure 
[5]). The rest of them belong to the quiet samples (i.e., flare-quiet ARs). The 
red diamonds denote the active samples, and the black dots represent the quiet 
samples in Figure [2J At the same time, the yearly mean values of the active 
and the quiet samples are calculated separately. Figure [3] shows these values. 
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Year Year 




1988 1993 1998 2003 2008 1988 1993 1998 2003 2008 

Year Year 



Figure 2. A<j>, Aip, |J Z |, |^c|i |«av|> Pfrcc; ^Ei an d ^Em of each AR sample during 1988-2008. 
Black dots represent the samples that did not produce flares with FI > 10.0 in the following 24 
h (flare-quiet samples). Red diamonds denote the samples that produced flares with FI > 10.0 
in the following 24 h (flare- productive samples). Two plots of monthly mean sunspot numbers 
during the same period are placed at the top of other panels for reference. 
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1988 1993 1998 2003 2008 1988 1993 1998 2003 2008 

Year Year 




1988 1993 1998 2003 2008 1988 1993 1998 2003 2008 

Year Year 

Figure 3. Yearly mean values of A</>, Ai/j, \ J X \, \hc\, |ct av |, Pfree , g(e> and fiEm of AR samples 
during 1988-2008. Dots represent the yearly mean values of the samples that did not produce 
flares with FI > 10.0 in the following 24 h (flare-quiet samples). Diamonds denote the yearly 
mean values of the samples that produced flares with FI > 10.0 in the following 24 h (flare-pro- 
ductive samples). The monthly mean sunspot numbers during the same period are overlapped 
(dashed line) for reference. 

The error bars indicate the standard deviations of the yearly mean values. The 
plot of monthly mean sunspot numbers (dashed line) is overlaid in each panel 
in Figure [3] for reference. 

As presented in Figure [2j A(j> and Aip have similar distributions. The linear 
correlation coefficient between them is 0.875. Neither of the two shear angles 
varies with the solar cycle. Most of their values are distributed in the range of 
[10°, 40°]. The points of Aip are more concentrated in this range. The mean 
value of ZY0 is 22.0° with the standard deviation of 6.6° during 1988-2008, 
while is 23.4° ± 7.6°. 
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Table 2. Linear correlation coefficients between eight parameters and the mean 
sunspot number (SSN), in yearly and monthly bases. 







Aip 


W 


\hc\ 


|<*av| 


Pfrcc 


''e 


(fern 


Mean SSN (yearly) 


0.379 


0.010 


0.093 


0.594 


-0.176 


0.666 


0.799 


0.814 


Mean SSN (monthly) 


0.242 


0.085 


0.054 


0.448 


-0.067 


0.474 


0.367 


0.530 



Combining with Figure[3j there is no obvious fluctuation of the parameter | J z \ 
with the phase of the solar cycle, in which most of AR samples have \J Z \ ranging 
in 1.0-3.0xl0 3 A km~ 2 . Some active samples at solar maximum have larger 
\J Z \. The samples with \J Z \ larger than 2.5xl0 3 A km -2 are likely connected 
with eruptive active regions. During the period from 2001 to 2003, which is the 
beginning of the declining phase of solar cycle 23, the values of \J Z \ are mostly 
lower than the general level of 2.0xl0 3 A km -2 . There are 91.7% of all the 
samples in 2001-2003 less than 2.0 xlO 3 A km" 2 . After the year 2003, the values 
of \J Z \ rise to the average level again. Figure [3] shows that the distributions of 
| J z \ is very similar in form to the two shear angles, A<f> and Aip. 

In the panel of |a av | , the values mostly concentrate in the range of 0-1.5 x 10 -5 
km -1 . There are no obvious changes in the parameter |a av | with the the solar 
cycle. The difference of |a av | between flare-productive and flare-quiet samples is 
insignificant except in the years 2004 and 2005. 

The parameters \h c \, pf rcG , and 4m follow nicely the variation of solar activity 
cycle. There are two bell-shaped parts in the corresponding panels during these 
two solar cycles. The overall level of these parameters rises gradually toward the 
solar maximum; while near the solar minimum they all fall down. Comparing 
with the monthly mean sunspot numbers, the troughs of 1989-1990 may be 
caused by the lack of the available samples or by some unknown uncertainties 
or evolution processes, but this result is consistent with |Bao and Zhang| ([1998 ) 
who have surveyed the evolution of the average current helicity in solar cycle 22. 
Table[Hlists the linear correlation coefficients between these eight parameters and 
the mean sunspot numbers in the yearly and monthly bases. These coefficients 
indicate that there is a much closer relationship between the mean sunspot 
number and \h c \, pf r0G , and 4m- 

Comparing the two parameters 4 and 4m representing complexity in mag- 
netic configuration, 4 of some ARs in the activity maximum periods is higher 
than that in the minimum periods, but the number of those ARs with higher 
4 are not many in the samples, and most of the samples arc distributed in a 
narrow band of ± 0.93 around 4 = 1.04. Combined with the average strength of 
magnetic fields, 4m reflects the magnetic complexity of an AR more accurately, 
while 4 is only related to the morphology but limited physics. Both 4 and 
4m change with the solar cycle, suggesting that there are more complex ARs in 
the activity maximum periods and less complex ARs in the activity minimum 
periods. 

Using magnetic synoptic charts from Michelson Dopplcr Imager onboard 
SOlar and Heliospheric Observatory (SOHO/MDI) and from National Solar 
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Observatory (NSO/Kitt Peak), IGuo et al] $ZUW\i found that the proportion 
of complex ARs (with c?e > 1) to all ARs decreases, although with a large 
fluctuation, in the declining phase of solar cycle 23, and the figure therein shows 
a small peak in the year 2005. Investigating the long-term evolution of magnetic 
hclicity in solar cycle 23, |Yang and Zhang] (|201ip found that the accumulation 
of magnetic helicity flux has the same trend in both hemispheres around 2005. 
|Zhang et aL|(|2010p and also |Tiwari, Ven katakrishnan, and Sankarasub ramanian| 
(|2009p found that ARs in the declining phases of solar cycles 22 and 23 do 
not follow the general hemispheric helicity sign rule. A similar result has been 
indicated by |Hao and Zhang| (|201ip . Figures [5] and [3] also show a small peak in 
2005 in each of the charts. Checking the monthly mean sunspot number and the 
monthly 10.7cm solar radio flux, there is no abnormal variation around 2005. 
There are only a little more flares in 2005 and in the end of 2004. Table [3] shows 
a comparison between the declining phases of solar cycles 22 and 23. We choose 
two years with roughly the same mean number of sunspots to compare the 
corresponding parameters, which means the solar activity level is similar in the 
two years of each pair. The exception is the mean sunspot number of 2002 and 
2003 to compare with the number of 1992. The quantities in Tablc[3]rcprcscnt the 
relative differences between the values of cycle 23 and cycle 22 ( X2:> ^ 22 1 where 
£23 and 222 represent the yearly mean values of the same parameter in the year 
of cycle 23 and cycle 22, respectively). The total activity level in the declining 
phase of cycle 23 is generally lower than that in cycle 22, and coincidentally the 
parameters in Table [3] follow the same tendency, that is, the parameter values 
of cycle 23 are lower than of cycle 22. However, all the parameters of the active 
samples have larger positive differences between 2005 and 1994. This shows that 
the nonpotentiality in flare-productive ARs is strong in 2005. A non-monotonic 
decline of cycle 23 might be regarded as one of the precursors to the long and 
deep minimum between solar cycles 23 and 24. 

4.2. Nonpotentiality Associated with Flares 

One of the generally accepted flare classifications is the SXR classification. From 
1975 to the present time, the GOES satellite has been recording the whole-Sun 
X-ray fluxes at 0.5-4 A (hard channel) and 1-8 A (soft channel) wavelength 
bands. The SXR flare classification (B, C, M, and X classes) utilizes the flux in 
the 1-8 A range based on the order of magnitude of the peak burst intensity. For 
analyzing the relationship between the nonpotentiality and solar flares, a time 
window r has been defined as the period forward in time (toward later times) 
from the observed vector magnetogram. Within a fixed time window r, the SXR 
flare index summed by weighting the flares of different classes is used to measure 
the flaring capability of an AR: 

FI = 100^/x + 10^/ M + ^/c, (11) 

T t t 

where 7x, Im, and Iq represent the indexes of X-, M-, and C-class SXR flare 
events, respectively ( jAntalova, 1996||Abramenko, 2005"] ) . Because the X-ray back- 
ground in the solar maximum is too high to detect B-class flares (c/., Fcldman, 
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Year AS 


QS TS Year AS QS TS 


02<< 


;03 vs. 92 -0.05 

04 vs. 93 0.08: 

05 vs. 94 0.40' 

06 vs. 95 0.16; 

07 vs. 96 0.29! 


' -0.398 -0.329 
i — U.ooo — U.zyo 
7 -0.183 -0.101 
> -0.170 -0.126 
i 0.031 0.070 


02&03 vs. 92 0.009 -0.269 -0.215 
U4 vs. yo U.z44 — U.2oy — U.zl4 

05 vs. 94 0.747 0.218 0.291 

06 vs. 95 0.372 0.068 0.108 

07 vs. 96 0.450 0.129 0.177 


02<< 

W\ 


;03 vs. 92 -0.35, 

04 vs. 93 0.09: 

05 vs. 94 0.59 

06 vs. 95 0.21: 

07 vs. 96 -0.03 


! -0.560 -0.509 

r — O.z/Z — U.Z4U 

_ 0.136 0.199 
! 0.086 0.102 
! -0.244 -0.212 


02&03 vs. 92 -0.603 -0.906 -0.753 

04 vs. \)6 — O.oOo — U.bUy — O.bOU 

05 vs. 94 0.074 -0.442 -0.370 

06 vs. 95 -0.208 -0.336 -0.320 

07 vs. 96 -0.493 -0.538 -0.531 


02& 

|c*av| 


;03 vs. 92 0.89 

04 vs. 93 1.571 

05 vs. 94 6.39! 

06 vs. 95 1.98( 

07 vs. 96 0.19: 


. -0.135 0.032 
i -0.236 -0.074 
i 1.325 2.028 
1 0.280 0.502 
! -0.042 -0.008 


02&03 vs. 92 -0.609 -1.001 -0.768 
04 vs. 93 -0.415 -0.696 -0.671 
p-— 05 vs. 94 0.687 -0.407 -0.256 

06 vs. 95 0.027 -0.467 -0.402 

07 vs. 96 -0.565 -0.502 -0.511 


02<< 


;03 vs. 92 0.16; 

04 vs. 93 -0.09 

05 vs. 94 0.45- 

06 vs. 95 -0.35' 

07 vs. 96 -0.26' 


> -0.282 -0.212 
1 -0.197 -0.189 
: -0.444 -0.320 
) -0.278 -0.288 
' -0.180 -0.193 


02&03 vs. 92 -0.105 -0.661 -0.521 
04 vs. 93 -0.425 -0.491 -0.485 
d Em 05 vs. 94 0.319 -0.637 -0.505 

06 vs. 95 -0.432 -0.455 -0.452 

07 vs. 96 -0.441 -0.384 -0.393 









flares are not included in the computation. FI = 10.0 (equivalent to a M1.0 
flare) is adopted to divide flare-productive ARs from flare-quiet ARs in Section 
14. H and also in the following parts. 

In Figure [3j the yearly mean values of each parameter for flare-productive 
ARs are higher than those for flare-quiet ones in the solar maximum periods, 
especially in the panels of \h c \, pf lcc , and <iEm- However, the twist factor |a av | 
shows an insignificant difference between the active samples and quiet samples, 
but shows more apparent difference with much higher uncertainty in 2005. For 
the same r, the yearly mean values of each parameter for the flare-productive 
ARs increase as the threshold of FI is increased, and the differences between 
flare-productive ARs and flare-quiet ARs are gradually enlarged accordingly. In 
reverse, for the same threshold of FI, the yearly mean values of each parameter 
for the flare-productive ARs decrease as the time window r is increased, and the 
differences between flare-productive and flare-quiet ARs are correspondingly a 
little reduced. Setting a series of r and different threshold of FI, the yearly-mean 
values (Figure[3]) will show basically similar distributions and evolution trends for 
each parameter. However, when the threshold of FI for active samples is raised 
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to a certain limit (e.g., M3.0 within 24 h; M5.0 within 48 h), the flare-productive 
samples vanish in some years of solar minimum periods, especially in the years 
1994-1996 and 2006-2008. It suggests that the flare-productive ARs are more 
likely to have relatively strong nonpotcntiality and great complexity, and those 
parameters characterizing nonpotentiality may be applied as indicators in the 
flare prediction. 

Furthermore, to study the flaring probability of an AR with certain properties 
of no npotentiality and complexity, the solar flare productivity is adopted ( Cui 
et aL 120061 120071 |Cui and Wang[ 120081 |Park, Chae, and Wang[ [ 2010p . which is 
defined as 

where N^(> X) is the total number of samples with the values of each parameter 
greater than its threshold X, and N&(> X) is the number of active samples with 
the values exceeding the same X. 

The flare productivities are plotted for each parameter in Figure |U To guar- 
antee the reliability of the statistics, the total number of samples is no less than 
50. In Figure |H the flare productivities are indicated by plus, cross, diamond, 
and asterisk symbols when the time window is set as 6, 12, 24, and 48 h from the 
time of the magnetogram. The Boltzmann sigmoid function is used for fitting 
the data. A(j> and Atjj have nearly the same tendency. \J Z \ is also similar to the 
above two, but rises faster at the value around 2.5xl0 3 A km -2 . \h c \ and |a av | 
have similar increasing tendency in this plot. The shape of the fitting curve of 
Pfrcc bears a resemblance to that of d^m- When the value of c?e is greater than 
1.0, there are almost no distinct differences in the flaring capability. These plots 
reveal different flare-productivity levels and their respective variation tendencies, 
as the result of the distributions and variations of the eight parameters in Figures 
[5] and [31 In Figure [H the fitting curves of the parameters with less differences 
between flare-productive ARs and flare-quiet ARs are relatively flatter than 
those with more differences. The largest flare productivity versus the largest X 
values in all the panels of Figure [3] may reflect their respective contributions to 
the flare prediction, and these specific parameters could be given corresponding 
weight in prediction models. 

4.3. Comparison among Nonpotentiality Parameters and Effective Distances 

The linear correlation coefficients for every two nonpotentiality and complexity 
parameters are listed in Table HI Being consistent with the results mentioned in 
Sections 14.11 and 14.21 c?Em has a relatively close relation with \h c \ and pf ICC ■ 4m 
has relatively higher positive correlation with nonpotentiality than <#e, while 
c?e have almost no relation with nonpotcntial parameters. The directly related 
parameters A<f> and A-0, c?e and c?E m , and \h c \ and \J Z \ arc highly correlated 
with each other as expected. The parameters \h c \ and pf ree also show very 
high correlation. Figure [5] displays scatter plots of these four pairs which show 
the highest correlation coefficients, 0.911, 0.875, 0.865, and 0.737. The flare- 
productive samples and flare-quiet samples arc plotted as red squares and black 
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Figure 4. Flare productivity for the parameters A<f>, Ai/>, \ Jz\, \hc\, |c*av|g Pfrccj an d efEm- 
The plus, cross, diamond, and asterisk symbols denote the forward time window of 6, 12, 
24, and 48 h from the time of the magnctogram, respectively. The solid lines are Boltzmann 
sigmoid curves fitted to the data. 
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Table 4. Linear correlation coefficients for eight nonpotcntiality and 
complexity parameters. 
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Figure 5. Correlations between A(p and Aip, \h c \ and Pfrccj 4 an d 4m, |^c| and |J 2 |. In 
each panel, the black crosses denote the flare-quiet samples, and the red squares denote the 
flare-productive samples; the solid line is the least-square quadratic polynomial fitting curve. 
'Linear CC stands for the linear correlation coefficient, and its value is obtained from all the 
data points of a panel. 
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crosses, respectively. The least-square quadratic polynomial fitting is made for 
all the data points of each panel. 

5. Conclusions 

By calculating eight parameters (A0, Aip, |J 2 |, \h c \, |a av |, pf xcc , d&, and dE m ) 
of nonpotcntiality and complexity for 2173 photospheric vector magnetograms 
in 1106 ARs associated with flares, we found the main results as follows: 

(1) On average, two mean magnetic shear angles Ac/> and Aip, mean absolute 
vertical current density \J Z \, absolute twist factor |a av |, and effective dis- 
tance cIe in ARs do not change significantly with the global solar activity 
level. However, it is more likely that these parameters show higher values 
in the solar maximum than in the solar minimum. 

(2) The mean absolute current helicity density \h c \, mean free magnetic energy 
density pf rco , and modified effective distance c?Em show high positive cor- 
relation with the mean sunspot number, and these parameters also have 
relatively close relationship with each other. The Pearson linear correlation 
coefficients of the above three with the yearly mean sunspot numbers are 
larger than 0.59. They can be used to characterize the solar activity level 
as well as the traditional sunspot number. 

(3) The nonpotcntiality and complexity parameters between flare-productive 
ARs and flare-quiet ones computed in this work arc useful to understand 
the evolution of flare-productive ARs and the relationship with the mag- 
netic activity levels of the cycles. These nonpotcntiality and complexity 
parameters may be synthetically applied as indicators to predict solar flares 
with some weight. 

(4) Due to the loss of the information of magnetic field strength in the parameter 
of effective distance cIe, the modified effective distance <iEm (including the 
strength of the magnetic field) turns out to be much better in indicating 
the magnetic activities of ARs. 

The accumulated nonpotential energy in the magnetic field of solar ARs pro- 
vides the energy of solar flares, while the trigger of flares is also related to the 
magnetic nonpotcntiality. The parameters discussed above will provide the basic 
information of nonpotcntiality of solar active regions in different aspects. The 
synoptic analysis of different nonpotential parameters may supply an objective 
basis on the relationship between the nonpotential field and solar flares in order 
to predict solar flares. 

Although long-term accumulation of the data from ground-based instruments 
is indispensable, the vector magnetograms taken by the Hinode satellite, Solar 
Dynamics Observatory (SDO), and other space-borne solar telescopes without 
being affected from atmospheric seeing have been giving more accurate data 
with better spatial resolution for further comparison studies. 
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